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Abstract 

A  linearized  three  component  MHD  theory  Is  given.   The  three 
components  are  two  oppositely  charged  species  and  a  neutral 
species.   Linearized  macroscopic  equations  are  presented  which 
in  the  limit  of  vanishing  magnetic  field  reduce  to  fluid 
dynamical  equations,  and  in  the  limit  of  vanishing  neutral 
component  give  the  Lundquist  equations.   The  coupling  of  the 
components  is  frictional.   In  addition  the  sound  speed  that 
enters  is  a  function  of  the  three  equilibrium  densities.   These 
equations  yield  a  dispersion  relation  which  is  tenth  order  in 
the  wave  speed   (to/k).   The  decay  frequencies  (i.e.,  -Im(a))) 
of  all  ten  modes  are  calculated  to  first  order  terms  in  the 
friction,  and  are  found  to  be  all  isotropic.   In  addition  all 
of  the  modes  are  stable. 


-  11  - 


NYO-9189 


Table  of  Contents 

Page 

Abstract  11 

Section 

I.   Introduction  and  Summary  of  Results  2 

II.   Analysis  5 

A.  Macroscopic  Equations  5 

B.  Dispersion  Relations  12 

References  22 


-  1  - 


NYO-9189 
Linearized  Three  Component  Magneto -Hydrodynamics 

I.   Introduction  and  Summary  of  Results 

This  analysis  Is  motivated  by  the  fact  that  most  experimental 
and  natural  (i.e.,  astronomical)  plasmas  are  rarely  fully 
Ionized,   In  this  paper  we  will  Initiate  a  study  of  the  magneto- 
hydrodynamlc,  or  low  frequency  waves  which  can  propagate  In  a 
three  component  plasma  In  a  steady  magnetic  fields   The  third 
component  Is  neutral. 

The  equations  of  motion  are  similar  to  those  stated  by 
Prank-Kamenetskll  [l],  Watnabe  [2j  and  Schliiter  [3]-   In  the 
first  two  works  the  extension  of  the  "cold",  two  component 
plasma  Investigations  [4,5],  to  the  three  component  case  was 
carried  out. 

The  present  work  Is  purely  macroscopic.   However  we  are 
guided  by  the  result  of  a  previous  microscopic  analysis  by  the 
author  [6]  In  which  all  of  the  long  wavelength  phenomena  of  a 
two  component  plasma  In  a  steady  magnetic  field  was  uncovered, 
and  formally  classified.   In  this  sense,  although  the  Included 
work  Is  strictly  adhoc,  It  Is  consistent  with  certain 
observations  partalnliig  to  the  transition  from  the  microscopic 
to  the  macroscopic  description  In  the  case  of  less  complicated 
two  component  fluid.   There  It  was  observed  that  the  six  MHD 
roots  [7,8,9]  (i.e.,  normal  modes  of  the  linearized  Lundqulst 
equations)  follow  from  the  more  general  dispersion  relation  In 
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the  limit  CO  =:  ak  +  .  .  .   (a  =  constant,   k  =  wave  number, 
CO  =  frequency).   Indeed,  It  Is  In  this  formal  limiting  process 
that  the  full  equation  of  motion  for  the  current  reduces  to 
the  simple   E  +  u  X  B   =  0  . 

0 

In  the  present  analysis  the  three  component  macroscopic 
equations  are  transformed  to  equations  which  describe  the  two 
component  charged  fluid,  (in  terms  of  charge  and  current 
densities)  Interacting  as  one  fluid  with  the  neutral  component. 
The  set  Involves      equations  resembling  the  Lundqulst  and 
neutral  fluid  [10]  equations.   The  frlctlonal  Interaction  appears 
as  forcing  terms  In  the  two  momentum  equations.   A  simple 
thermodynamic  treatment  yields  a  form  for  the  three  component 
sound  speed.   It  appears  as  a  function  of  Individual  equilibrium 
mass  densities  and  the  equilibrium  temperature.   The  model  Is 
such  that  for  vanishing  magnetic  field,  there  Is  only  one  sound 
speed,  which  is  believed  to  be  Indicative  of  the  true  fluid 
dynamical  (i.e.  macroscopic)  region.   The  equation  of  motion  for 
the  current  Is  kept  In  the  above  Infinite  conductivity  form, 
which  Is  consistent  with  the  af ore-mentioned  classification 
process . 

A  normal  mode  analysis  Indicates  that  the  dispersion  relation 
that  this  system  of  equations  yields  Is  tenth  order  in  co  .   In 
the  limit  of  vanishing  coupling  these  ten  roots  are  the  six 
MHD  roots,  the  two  sound  modes,  and  a  degenerate  co   factor 
which  is  the  limiting  form  that  the  "particle  paths"  dispersion 
relation  takes  in  the  present  linearization.   To  first  order 
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terms  In  the  friction  frequency  v  ,  the  decay  frequency  (i.e., 

-lm(co))  of  the  three  perturbed  modes  are  all:   vp^/2r   where 

p^   is  the  equilibrium  mass  density  of  the  neutral  component 

and  r   is  the  total  equilibrium  mass  density.   The  isotropy 

of  these  MHD  decay  frequencies  is  a  reflexion  of  the  purely 

frlctional  coupling  of  the  charged  and  neutral  species;  and  the 

null  effect  of  the  steady  magnetic  field  in  the  neutral  component. 

The  decay  frequency  of  the  perturbed  sound  mode  is 

vp  /2r    ,   where   p    is  the  mass  density  of  the  combined 
^o/   o  "^o 

charged  components.   Finally,  one  of  the  two  "particle  path" 
transverse  modes  remains  as  a  degenerate  od  factor,  while  the 
second  one  decays  four  times  faster  than  the  above  sound  mode. 
The  degenerate  root  indicates  that  to  terms  of   0(v)  the 
related  perturbation  does  not  decay. 

All  ten  modes  are  found  to  be  stable. 
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II.   Analysis . 

A.   Macroscopic  Equations. 

The  system  we  will  consider  is  a  three  component  fluid. 

Two  of  the  components  are  oppositely  charged  and  the  third 

component  is  neutral.   We  wish  to  examine  the  small  amplitude 

perturbation  of  the  system  about  a  certain  equilibrium  state, 

restricting  the  investigation  to  the  magneto-hydrodynamical 

(MHD)  region  of  phenomena,  wherein  the  density  of  the  charged 

components  is  large  enough  so  as  to  insure  fluid  dynamical 

reactions.   More  formally  this  constraint  appears  as  a  limit  on 

the  disturbances  which  will  propagate;  their  frequencies  oo 

being  such  that  to  «  cd   ,   where  co   is  the  characteristic 
^  o  o 

plasma  frequency  [ll]. 

The  neutral  component  is  similarly  constrained  in  as  much 
as  the  related  phenomena  are  governed  by  the  fluid  dynamical 
equations,  which  are  a  valid  description  for   A  »  A   ,  where 
A   is  the  mean  free  path  of  the  neutrals,  and   A   is  the 
wavelength  of  the  propagating  disturbance. 

As  stated  above,  the  analysis  describes  the  two  component 
charged  species  interacting  as  a  single  fluid  with  the  neutral 
species.   This  formulation  is  brought  about  through  the 
following  transformation  of  variables.   The  notation  is  such 
that   [p-,  ,  p  ,  p^]   are  the  three  mass  densities  where  the 
first  two  indices  relate  to  the  charged  components.   The  bar 
signifies  the  "full"  state  variable,  as  opposed,  say,  to  the 
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equilibrium  value,  eg. ,  p    and  the  perturbation  value   p  . 
The  three  are  related  through  the  equation   p  =  p   +  p   , 
which  in  turn  defines  the  perturbation   p  . 

In  similar  manner  the  macroscopic  velocities  are 
[u, ,  u  ,  u^]  .   The  first  two  variables  of  both  sets  are 
transformed  to  the  one  fluid  variable  according  to: 


(1) 


^1'  ^2   ^  ^'  ^ 


where 


Pi'  Pg  ^  P- 


-  _  Pl^l  ^  P2^2 
Pi  +  P2 


Pi  - 


PqU 


2  2 


J  =  e  —  Ut  -  e 

m.,   1      mp 


P  =  Pi  +  P2 


_  ^  Pi    ^  ^ 
^  m-j^    ^  m^ 


The  neutral  variables   (p^^u^)   are  kept  intact,  as  are  the 
temperatures   l^'^-i  j  '^p^  T  ]l. 

The  equilibrium  configuration  is  given  by  : 


J  =  u   =0 


o     o 
Pi     Po 


m-, 


m. 


o    o 


o    o 


•>      9      ^   9         '   P3  "  P3 


(3) 


o    0,0 
r  =  p   +  p 


3 


u 


3 


0   ;    T°  =  T°  =  T°  =  T°   ;    tE,B]^  =  [0,B^] 
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Note  that  the  new  variable  r  Is  the  equilibrium  total 
mass  density.  As  mentioned  above,  the  perturbation  is  given 
by, 

(4)  J-J°+J;       u=u°+u,... 

The  conservation  equations  that  these  (full)  variables 
satisfy  are  [6]  : 

The  continuity  equations 


a    - 


P  +  \7  •  u  p  =  0 


(5)  .         1^  Q  +  V  •  J  =  0 


^  P3  +  y  •  p-ju^  =  0 


The   momentum   equat.ions 

_;^__  _  ______       A 

P^ft   +   u-  ^-O^  +   7-  Pi  j    -   Q(E  +  u   KB)    -    J  X  B  =  M^ 


a 


■ij 


(ft  +  v-u)j  +  jQ%^  +  {1  +  m)-  V  u  +  V-^ij 

f  c  \  2  ^  e(m-,  -m_  ) 

^^'  -    p   ^—    (E   +    u    kB) i— ^    (EQ    -    J   X  B) 


m-j^m^    ^  m-j^m^ 


P3(|t  +  Uj-V)u3   .V^''  =  A3 
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The  tensor  p.  .   is  the  pressure  term,  and  the  variable 
R.  .   represents  the   "J-Flux"  tensor  [6]  .   The  terms  /\  . 
are  functionals  of  the  collision  Integrals.   By  conservation 
of  momentum  in  a  collision  /  X-,  +  /  \  ^  ==  0  . 

The  above  equations  are  standard  moment  equations.   They 
are  written  down  here  to  indicate  the  effect  that  the  formal 
change  of  variables  has  on  the  equations  involved.   The 
succeslve  higher  order  moment  equations  (e.g.,  pressure,  heat 
flow,  . . . . )  would  be  of  the  same  form,  the  first  two  in  each 
case  relating  to  the  charged  components  and  the  third  to  the 
neutral  component.   It  is  clear  that  the  formal  change  of 
variables  yields  equations  which  describe  the  charge  components 
as  a  coupled  fluid  and  the  neutral  component  as  another 
interacting  fluid.   The  suggestion  is  something  less  than 
fortuitous,  as  the  terms  which  enter  into  the  charge  component 
equations  are  physically  relevant  (e.g.,  J,  Q).   These  are  the 
variables  which  interact  with  the  fields  in  recognizable  manner, 
which  in  turn  is  the  major  requirement  that  a  variable  must 
satisfy  in  an  ad-hoc  theory. 

If  the  above  equations  are  linearized  about  the  equilibrium 
state  as  given  by  equation  (3),  and  the  fluid  dynamical 
constraint  is  imposed  [6],  then  the  lowest  order  terms  (in  the 
wave  number  k)  appear  as  (including  Maxwell's  equations)  : 


P    +    PQ   \/"^    =     0 


P3  +  P3  V '^3  =  0 

p^u  +  a  V  p  =  J  x,  B^  -   ^  -^   (u  -  u  ) 

o 
(7)  o 

p  u  +  a  x/  p^  =  — ^r —  (u  -  u  ) 


E  +  uXB=0;         \/-B=0 
VXB=h^J      ;         ^XE=-B 

The  collision  forms  /\    and  /  \^   are  replaced  by  the 
momentum  exchange  terms  :   -   and   +   (p  p^/r  )v(u-u^)   , 
respectively.   The  single  sound  speed   a  ,  that  enters,  is 
consistent  with  simple  thermodynamic  considerations.   First 
we  recall  the  adiabatic  formula   P  =  Ar'  ,  where   P   is  the 
total  pressure,   2p.  ,   r   is  the  total  mass  density  2p .  , 
A   is  a  constant,  and  y   =    (d+2)/d  =  5/3   for   d  =  3,   i.e., 
three  degrees  of  freedom.   Each  component  obeys  an  equation 
of  state   p.  =  p.R.T   (the  thermodynamic  equilibrium  is  such 
that   T.  =  T  )  .   It  follows  that   P  =  r  R  T   where, 

(8)  R .  ^i!i_!j2!^^_p3!i 

For  the  sound  speed   a  ,  we  have  : 

(9)  -'  '^'y^-y^^o'^^^o■ 
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In  the  specific  configuration  under  considerations: 

(10)  R  =  K  [2n°  +  n°]  /[Mn°  +  m^n°]  , 

where  we  have  written  M  for  m.,  +  m   ,  and  have  recalled  that 
the  equilibrium  configuration  Is  charge  neutral,  so  that  the 
equilibrium  number  densities   n-,   ,  satisfy  n-,  =  n^  .   k   Is 
Boltzmann's  constant. 

Before  examining  the  linearized  system  (7)  for  plane  wave 
solutions  let  us  first  consider  two  limiting  cases. 


)    B^^O 


In  this  Instance  the  system  (7)  reduces  to 


r_U  +  a^^r  =  0 


o 


oV 


r  +  r  W'U  =  0 


It    t-n    -V    "-"5^=° 


Capital   U  represents  the  mass  averaged  velocity 

o    ,0    ,0 
p  u  +  p  u^  +  p^u^ 

-^  -^    ^ i— ^  ,   while   r  represents  the  total  mass 


r 
o 


density  p.,  +  p„  +  p 


3 
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In  the  limit   B  >  0  ,   the  coupling  to  Maxwell's 

equations,   J  .  B    ,  vanishes.   Furthermore  since   B  =  0  , 
^        '        o  o 

Ohm  law  gives   E  =  0  ,  which  In  turn  Implies  that   B   Is 
constant  in  time.   The  only  remaining  relevant  field 
equation  is  \/  X  B  =  J  ,   which  is  decoupled  from  the  set 

(11). 

The  set  (ll)  yields  the  dispersion  relation: 

(12)      cn^(cD^-a^k^)(ao+lv  )^[a)(aj+iv  )  -  a'^k^]  =  0   . 

The  first  and  third  factors  give  the  eigenvalues  of  the 
transverse  components  of   U  ,   and   u-u^  ,  respectively.   For 
U  say,  these  components  may  be  written   k  x  U  ,   and 
k  ,>(  (U  Y  k)  .   The  second  factor  gives  eigenvalues  for  the 
longitudinal  component   k-U  ,   of  U  ;  and  for  the  total  mass 
density  r  .   The  last  factor  relates  to  the  longitudinal 
component  of  the  difference  flow   (u-u^)  ,  and  the  difference 
mass  density   [(p/p  )  -  (p^/P^)]-   It  should  be  noted  that  in 

the  space  homogeneous  limit   k  )  0  ,  this  last  factor  yields 

the  two  roots  co  =  0  and   ito  =  v  jwhich  relate  respectively 
first  to  a  state  wherein  the  perturbation  difference  in  mass 
densities  remains  constant,  and  second  to  state  wherein  the 
two  flow  velocities  relax  to  each  other.   The  decay  time  for 
this  latter  phenomena  is   v~   .   More  generally  the  third 
factor  relates  to  similar  relation  of  the  transverse  vibration 
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(b)    P3  — >  0 


In  this  second  limit  the  system  (7)  reduces  to 


p  ^  po\7 


u  =  0 


(13)  p^ix  +  a^  S^p   =    J  >^  B^ 


E  +  u  X  B   =  0  , 
o      ' 


together  with  Maxwell's  equations.   This  aggregate  of  equations 
is  recognized  to  be  the  standard  linearized  Lundquist  set,  which 
yield  the  so  called  Past,  Slow  and  Alfven  modes  [7-9]- 


B.    Dispersion  Relations 

We  will  now  examine  the  full  system  (7)  for  normal  mode 

solutions.   First  we  note  that  the  level  of  description  of  the 

Included  analysis  is  given  by  the  dimensionality  of  the  state 
vector   I   , 


(1^) 


)  E,B,u, J,p,u^,p^^ 


The  independent  components  of  \~^     are  as  follows:   three 
components  of   u  ,  two  components  of   B  ,  three  components  of 
u   ,  together  with  p   and   p   .   These  components  Individually 
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satisfy  first  order  time  differential  equations.   The  time 
dependance  of  the  remaining  variables  are  determined  in  terms 
of  these  independent  variables. 

The  plane  wave  transform   |  '    of  the  state  vector   f 
Is  defined  as: 

(15)  I   (x,t)  =  I  '  (co,k)e  ^ 

In  terms  of  this  state  vector   |    ,  the  system  appears  as: 

(16)  A,.(|^  '  1^  '  V  ;   p^,p;,a2,B^,c,sjnj  =  0 

The  transform  of  the  above  equation  appears  as  (writing 
only  the  dependence  of  the  first  three  parameters): 


(17)  ^^j^^'  ^'  ''^^j*  "  ° 


which  yields  the  dispersion  relation 


(18)  D  =  lA^j(a3,k,  v)|  =  0 


The   function     D      is   a   tenth   order   polynomial   in     (Jo    .      If 

loose      B     =    (0,0, B    )      ai 
o        ^    '    '    o  ^ 

as  (after  some  elimination) 


we  choose   B  =  (0,0, B  )   and   k  =  (k  ,0,k  ),   then  D  appears 
00  X     z 
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(19) 
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where , 


(20)      V  =  (p^/r^)v    ;     v  =  (p  /r^)v 


We  will  now  proceed  to  find  the  roots  of  the 
dispersion  relation  (l8)  to  terms  of   0(v).   To  these  ends 
D   is  expanded  in  a  Taylor  series  about   v  =  0   : 


I  2     I! 

(21)      D(v)  =  D(0)  +  vD  (0)  +  -^  D  (0)  + 


2 


where   D  appears  formally  as  : 


(22)      D=^  (-)%e/2e/3e3  ••'  ^6e^ 
( s-j^,  e^,  e^,  •  .  •  Sg  ) 


where   p   is  the  order  of  the  permutation   P(e-,,e  ,  ...)•   It 
follows  that  the  derivative  of  D   is  given  by  : 


(23)      d'  =  21  (-)^d|g  d^g  d    ...  +  Z_  (-)%-  d..  d 


le^  2e^   ^e    '  '  '         ^ —      le-^^  2e^   ^e^ 


+  SI  (-)^d^e  '^2e  ^^   •••  +  •• 


.e-^  2e2  3e^ 

The  n   sum  is  the  original  determinant  with  the  n   row 
replaced  by  the  derivative  of  that  row.   If  these  determinants 
are  evaluated  and  then   v   is  set  equal  to  zero  It  Is  easily 
verified  that  the  only  terms  that  contribute  to  D  (o)   are 
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the  diagonal   3  \  3   subdeterminants .   It  follows  that  for 
purposes  of  evaluating  D( v )   to  terms  of   0(v),  we  may  set 


2h)  D  =  k^°D  D, 

o  1 


where 


(25) 


26) 


D  =  [v^-v^{a.^+k^)+a.^k^-v   v(a2+A2-2v^  )  ]  ( A^-v^-v  v  ) 

O  Z    V  Zi        V 


D,  =  (v  +  V  )^[a^-  v(v  +  V  )] 

1     ^        V  V 


The  new  variables  are  defined  as: 


v„  =  ^  ;   V  ^  1  E  ;   v^  =  i^  ;   k^  =  R  cos 

dv     p         ov     p^ 
V    "^o         V  _  "^3 


Sv     r    '     'bv  r^ 

V     o  V     o 
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The  Taylor  series  of  D(v  )   about   v  =0  now  appear  as: 


k-10D(vJ  =  D(0)  +  (D^  — ^  — ^  +  D^  — i  — ^)  V   + 


V     V 


V     V 


Vg  =  0 


O     V 
V     V 


o 


-^  V [ 2a^A2+A^A2_2v^A2-2v2A^-2v^a^+3v^ ] 
o 


V   =  0 

V 


1    V 


bv      bv 

V     V 


=  ^  V  [2a2-  3v2] 


V  =  0 

V 


The  unperturbed  sub-determinants   D  (o)   and   D  (o)  are 


D^(0)  =  [v^-v^{a^+A^)   +   a^A^]  (A^-v^) 


(28) 


D-^(O)  =  (a2-v2)v2    , 


where  we  recall 


(29)      D(0)  =  D^(0)  D^(0) 


The  roots  of  D  (o)  are  the  six  modes  of  the  Lundqulst  [8] 
relation,  viz  : 

(30)       D^(0)  =  (a2-v2)  (v^-A^)  (v2-a2)    . 
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These  three  sets  of  roots  are,  respectively; 


(l)  The  Alfven  wave 


(31) 


v^  =  A^  =  A^cos^ 
z 


2)   The  "Fast"  wave 


(32) 


2    „2    f  a^+A^ 


v"  =  A' 


+ 


+ 


N 


2   2  2 
— -  ^  ^z 


(3)   The  "Slow"  wave  : 


(33) 


v2  =  a2  = 


a^+A^ 


a^+A^      2,2 
2      "  ^  ^z 


The  roots  of  D-|(0)   are  the  sound  mode  and  the  degenerate 
transverse  mode. 

It  remains  to  uncover  the  extension  of  these  roots  to  terms 
of   0(v  )  .   All  of  these  are  of  the  form  : 


:34) 


V.  =  V.  +  a.  V 

1      1      IV 


where   V   is  the  ten  component  vector  : 


35) 


V  =  (±A^,  +A^,  +A_,  +a,  ±0) 


Substitution  of  (3^)  into  the  dispersion  relation  (24)  and 
retaining  terms  to   0(v  )   give  the  desired  results. 
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First  we  note  that  the  D  (o)  D^fo)   factor  Involves  terms 

o     1 

?   p 
of  the  form  v.-V.   ,  so  that  the  difference  between  the  two 

(+)        J.  '    f+) 

substitutions        v.  '  =  IV.  +  a.  ^v   into  this  factor  is 

a  change  in  sign.   The  remaining  first  derivatives  terms 

accept  only  the  lower  ±V.  ,   part  of  v.  ^   since  a     part  would 

give  an   0(v  )   term.   In  addition  both  of  the  higher  order 

derivatives  terms  are  of  the  form  vf(v  )  ,  so  that  the  difference 

between  the  related  two  mirror  substitutions   tV.  ,      also 

involves  only  a  change  in  sign.   It  follows  that  the  same  equation 

results  for  a.      and  a~    ,    so  that  it  is  consistent  to  write  : 
1        1 

(36)  v~  =  +V.  +  a.v 

^■^  ^  1     1  V 


This  observation  together  with  the  observation  that 
a.  <  0  ,   is  merely  a  consequence  of  the  stability  of  the 
equilibrium  configuration. 

The  said  substitutions  yield  the  following  results: 


I  :    The  MHD  Roots 


(37) 


(a)   Alfven  Modes 


v.,  =  +A  +  a^v 
1    -  z     1  V 


°1  =  -  P3/2^o 
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(b)   Fast  and  Slow  Modes 


V 


+ 


2,1 


^^-.  +  <\ 


(38) 


V   +  =  ±A   +  o^v 


2  V 


where 


(39)    2a 


(±)  _  P3 


o  r„  2 ,  ^  2 .  „  2 


A^(2a^+A^ 


2A|(A^+A^+a^)  +  3A+ 


o   L  (A^  -  a2)  (a2  -  A?) 


If  the  values  for  A   ,  A,  as  given  by  (31-33)  are 
substituted  Into  the  above^  one  finds  that 


40) 


^o  =-pS/2r. 


3/^0 


The  decay  times  of  the  three  MHD  waves  are  all  equal  and 
Isotropic . 

II     The  Neutral  Gas  Roots 

The  extended  sound  modes  are  written 


(^1) 


vj  =  ±a  +  a^v^  , 


which  when  substituted  Into  equation  (24)  yields 


This  observation  is  due  to  N.  Ainbender,  of  the  CIMS 
computational  group. 


-  20 


h2)  o     =    -    p  /2r 


3     ^o/"  o 


One  of  the  remaining  transverse  modes  remains  zero  to 
0(v  )  (i.e.  the  perturbation  remains  In  a  steady  state)  while 
the  remaining  mode  gives : 


(43)      v.,  =  -2v   p  /r    , 


a  very  fast  decaying  root. 
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